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Abstract 


Atomic  force  microscopy  performed  on  single  crystals  of  three  different 
polymorphs  of  bovine  insulin  revealed  molecularly  smooth  (001)  layers  separated  by  steps 
whose  heights  reflect  the  dimensions  of  a  single  insulin  hexamer.  Whereas  contact  mode 
imaging  caused  etching  that  prevented  molecular-scale  resolution,  tapping  mode  imaging  in 
solution  provided  molecular-scale  contrast  that  enabled  determination  of  lattice  parameters 
and  polymorph  identification  while  simultaneously  enabling  real-time  examination  of 
growth  modes  and  assessment  of  crystal  quality.  Crystallization  proceeds  layer-by-layer  in 
which  the  protein  molecules  assemble  homoepitaxially  with  nearly  perfect  orientational  and 
translational  commensurism.  Tapping  mode  imaging  also  revealed  insulin  aggregates 
attached  to  the  (001)  faces,  their  incorporation  into  growing  terraces,  and  their  role  in  defect 
formation.  These  observations  demonstrate  that  tapping  mode  imaging  is  ideal  for  real-time 
in  situ  investigation  of  the  crystallization  of  soft  protein  crystals  of  relatively  small  proteins 
such  as  insulin,  which  cannot  withstand  the  lateral  shear  forces  exerted  by  the  scanning 
probe  in  conventional  imaging  modes. 
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Introduction 


Correlation  of  protein  structure  with  biological  activity  relies  heavily  on  the 

determination  of  the  protein's  three-dimensional  structure,  typically  by  X-ray  diffraction 

which  can  reveal  protein-receptor  binding  site  conformations  and  interactions  that  are 

relevant  to  the  rational  design  of  pharmaceuticals  (De  Vos  et  al.,  1992;  GUnther  et  al., 

1990;  Lambert  et  al.  1989).  However,  one  of  the  primary  obstacles  to  protein 

characterization  is  the  growth  of  crystals  suitable  for  X-ray  diffraction.  This  difficulty 

stems  from  numerous  factors,  including  convective  currents  which  lead  to  poor  crystal 

quality  (Koszelak  et  al.,  1995;  Day  and  McPherson,  1992;  DeLucas  et  al.,  1989;  DeLucas 

et  al.,  1991;  DeLucas  et  al.,  1986;  Erdmann  et  al.,  1989),  and  aggregation  processes  which 

can  lead  to  non-crystalline  solids  or  crystals  with  high  defect  density  and  irreproducible 

solvent  contents  (Matthews,  1968).  Crystal  purity,  size,  morphology  and  polymorph 

identity  also  are  important  issues  in  the  manufacturing,  ease  of  delivery  and  bioavailabUity 

of  crystals  of  therapeutic  proteins  (Byrn,  1982;  Carstensen,  1977),  requiring  careful 

optimization  of  crystal  growth  conditions  such  as  temperature,  pH,  concentration,  and 

additive  levels.  The  need  to  improve  protein  crystallization  has  prompted  experimental 

investigations  (Georgalis  et  al.,  1993;  Sazaki  et  al.,  1993;  Thibault  et  al.,  1992;  Auersch  et 

al.,  1991;  Kam  et  al.,  1978;  Durbin  and  Feher,  1990)  and  modeling  of  the  nucleation  and 

growth  of  protein  crystals  aimed  toward  understanding  critical  nuclei  dimensions,  relative 

rates  of  nucleation  and  growth,  macromolecular  aggregation,  and  the  attachment  of  single 

protein  molecules  and  their  aggregates  to  growing  crystals  (Oosawa  and  Kasai,  1962; 

Noever,  1995a,  1995b).  These  issues  become  particularly  important  in  therapeutically 

important  proteins  such  as  insulin  (Mr  =  5.8  kDa),  the  hormone  responsible  for  regulating 

glucose  metabolism.  The  absence  of  this  hormone  leads  to  diabetes,  a  disease  which 

affects  approximately  1  out  of  every  10  persons  in  North  America.  First  crystallized  in 

1925  (Abel  and  Geiling,  1925;  Abel  et  al.,  1925;  Abel,  1926;  Abel  et  al.,  1927),  insuUn 

single  crystals  consist  of  ordered  arrangements  of  hexamers  assembled  from  three  insulin 

Yip  and  Ward  3 

Submitted  to  Biophysical  Journal 
March  27,  1996 


dimers.  The  bioavailability  and  efficacy  of  microcrystalline  insulin  medications  depend 
upon  the  crystal  dissolution  rate,  which  can  be  influenced  by  crystal  polymorphism, 
morphology  and  size,  and  the  stability  of  the  hexamer  toward  dissociation  to  active 
monomer  form  (Hollenberg,  1990;  Gammeltoft,  1988;  Berson  and  Yalow,  1966;  Smith  et 
al.,  1984).  Optimization  of  insulin  crystallization  would  be  aided  significantly  by  direct, 
real-time  visualization  of  nucleation  and  growth,  particularly  with  regard  to  observing 
critical  nucleation  events,  examining  the  role  of  aggregates  and  defects,  identifying 
polymorphs,  assigning  actively  growing  crystal  faces,  and  determining  crystal  quality 
under  process  conditions.  These  needs  prompted  us  to  examine  insulin  nucleation  and 
growth  by  real-time  in  situ  atomic  force  microscopy  (AFM),  which  has  been  used 
successfully  for  examining  the  nucleation  and  growth  of  a  variety  of  soft  organic  (Hillier 
and  Ward,  1994;  Hillier  et  al.,  1994;  Carter  et  al„  1994)  and  protein  crystals  (Land  et  al., 
1995;  Malkin  et  al.,  1995;  Konnert  et  al.,  1994;  Durbin  et  al.,  1993;  Durbin  and  Carlson, 
1992). 


Materials  and  Methods 

Bovine  insulin  single  crystals  were  grown  using  protocols  provided  by  Dr.  G. 

David  Smith  (Hauptman-Woodward  Research  Institute.  Personal  communication).  The  R6 

rhombohedral  crystallization  protocol  follows:  Bovine  insulin  (Sigma,  St.  Louis,  MO)  was 

dissolved  in  0.02  N  HCl,  0.15  M  zinc  acetate  and  0.20  M  sodium  citrate  buffer.  To 

ensure  complete  dissolution,  0.75  N  NaOH  was  added  followed  by  5%  phenol  and  1.0  M 

NaCl.  The  solution  pH  was  adjusted  to  8.5  with  0.75  N  HCl.  After  standing  at 

approximately  38“C  for  3  days  and  subsequent  cooling  to  ambient  temperature,  small  clear 

rhombohedral  crystals  formed  after  several  hours. 

Atomic  force  microscopy  images  were  obtained  on  a  Digital  Instruments 

NanoScope  III  MultiMode  scanning  probe  microscope.  Contact  mode  imaging  was 
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performed  with  Si3N4  cantilevers  having  a  measured  spring  constant  of  0.074  N  m-1,  as 
determined  by  an  end-mass  resonance  technique  (Cleveland  et  al.,  1993).  Tapping  mode 
imaging  was  performed  with  identical  cantilevers  but  with  oxide-sharpened  tips  (Digital 
Instruments  Inc.  Santa  Barbara,  CA).  The  reported  AFM  images  were  taken  using  the  E 
scanning  head  which  has  a  maximum  lateral  scan  window  of  14.6  x  14.6  |im.  For 
imaging  in  air,  crystals  were  transferred  by  pipette  onto  either  a  polished  AFM  mount 
coated  with  a  small  amount  of  spray  adhesive  or  a  small  tab  of  freshly  cleaved  mica  affixed 
to  the  AFM  mount.  The  mother  liquor  was  then  quickly  wicked  away  and  the  AFM  tip 
positioned  over  a  suitable  crystal.  The  regions  imaged  were  free  of  macroscopic  surface 
defects,  as  determined  by  optical  microscopy.  Imaging  in  solution  was  carried  out  in  a 
combination  contact/tapping  mode  liquid  cell  (Digital  Instruments  Inc.  Santa  Barbara,  CA). 
The  insulin  crystals  in  their  mother  liquor  were  introduced  into  a  well  which  was  defined 
by  an  0-ring  placed  on  an  AFM  mount  coated  with  a  thin  layer  of  vacuum  grease.  The 
AFM  liquid  cell  was  carefully  placed  over  the  0-ring  and  loose  insulin  crystals  were 
removed  from  the  field  of  view  by  gentle  infusion  of  mother  liquor. 

Tapping  mode  imaging  in  solution  was  performed  at  cantilever  drive  frequencies 

ranging  from  6  kHz  to  32  kHz.  Viscous  coupling  between  the  cantilever  and  the  fluid 

results  in  a  number  of  broad  resonant  peaks,  and  therefore,  the  choice  of  an  appropriate 

drive  frequency  cannot  be  determined  a  priori.  The  cantilever  drive  and  setpoint  voltages, 

which  determine  the  cantilever  oscillation  magnitude  to  be  controlled  by  the  feedback  circuit 

during  imaging,  have  a  dramatic  effect  on  image  quality.  These  parameters  are  chosen  in 

each  experiment  such  that  image  quality  is  optimized  with  respect  to  resolution,  contrast 

and  sample  integrity.  Increasing  in  the  setpoint  voltage  in  tapping  mode  corresponds  to 

increasing  the  average  sample-tip  separation.  Small  increases  (typically  ca.  0.03  V)  in  the 

setpoint  voltage  from  an  optimized  value  during  insulin  imaging  results  in  deterioration  of 

image  quality  due  to  tip  disengagement  and  diminished  resolution.  Based  on  the 

aforementioned  cantilever  spring  constant  and  the  cantilever-photodiode  sensitivity,  this 
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suggests  that  the  imaging  force  is  approximately  100  pN  during  imaging,  in  qualitative 
agreement  with  previous  reports  (Radmacher  et  al,  1994).  The  optimum  cantilever  drive 
and  setpoint  voltages  will  depend  strongly  on  the  cantilever  spring  constant,  drive 
frequency,  and  coupling  between  the  tip  and  surrounding  fluid.  It  should  also  be  noted  that 
a  similar  decrease  in  setpoint  voltage  from  an  optimized  value  results  in  mechanical  etching 
of  the  sample.  Image  analyses  were  performed  on  low-pass  filtered  AFM  data  using 
Nanoscope  III  version  4.21b2  (Digital  Instruments  Inc.  Santa  Barbara,  CA.)  and  NIH 
Image  version  1.59  (National  Institutes  of  Health  -  available  by  anonymous  FTP  from 
zippy.nimh.nih.gov)  software  packages. 


Results 

AFM  images  of  the  (001)  faces  of  the  R6  rhombohedral,  T6  rhombohedral  and  R6 

monoclinic  polymorphs  of  bovine  insulin,  acquired  in  conventional  contact  mode  directly  in 

mother  liquor,  revealed  molecularly  smooth  terraces  with  micron  scale  lateral  dimensions 

(Figure  1).  The  heights  of  the  steps  separating  the  terraces  were  43  A  and  50  A  for  the 

rhombohedral  and  monoclinic  forms  (The  AFM  data  for  the  T6  rhombohedral  form  are 

indistinguishable  from  that  of  the  R6  rhombohedral  form,  and  therefore  are  not  reported  in 

detail  here).  These  heights  were  consistent  with  the  c-axis  dimensions  of  these  forms 

(Baker  et  al.,  1988;  Derewenda  et  al.,  1989;  Smith,  1995;  Smith  and  Dodson,  1992; 

Ciszak  and  Smith,  1994;  Kaarsholm  et  al.,  1989;  R6  rhombohedral :  space  group  R3,  a  = 

b  =  79.92  A,  c  =  40.39  A  as  indexed  on  hexagonal  axes  ;  R6  monoclinic  :  P2i  space 

group  with  lattice  parameters  a  =  61.23  A,  £>  =  61.65  A,  c  =  48.05  A,  P  =  1 10.50  ,  Smith, 

G.D.  Personal  communication).  These  values  are  in  good  agreement  with  the  estimated 

thickness  of  an  isolated  insulin  hexaraer,  which  has  been  described  as  an  oblate  spheroid 

with  a  diameter  of  approximately  50  A  and  a  height  of  35  A  (Blundell  et  al.,  1972).  In 

contrast,  AFM  images  acquired  in  air  using  contact  and  tapping  mode  revealed  terraces 

separated  by  steps  with  heights  ranging  from  23  A  -  29  A.  The  smaller  heights  in  air  are 
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attributed  to  dehydration  of  the  (001)  protein  layers  upon  removal  of  the  crystals  from 
solution.  The  observed  decrease  in  the  thickness  of  the  (001)  layers  upon  dehydration 
corresponds  to  a  «  40  %  volume  change  when  shrinkage  is  assumed  to  be  along  the  c-axis 
and  is  in  good  agreement  with  the  estimated  30%  -  35%  solvent  content  for  these 
polymorphs  (G.D.  Smith.  Personal  communication.;  Gursky  et  al.,  1992;  Badger  and 
Casper,  1991;  Dodson  et  al.,  1978).  The  hydrated  and  dehydrated  surfaces  also  differ  with 
respect  to  the  integrity  and  roughness  of  the  terraces,  which  fracture  upon  dehydration,  and 
the  viscoelastic  response  of  the  surface  protein  layers  to  applied  tip  forces  (Yip  and  Ward, 
in  preparation). 


[Figure  1] 

Despite  small  imaging  forces  estimated  at  30  nN  and  4  nN  in  air  and  solution, 
respectively,  all  attempts  to  obtain  molecular-scale  contrast  of  the  (001)  layers  using 
contact  mode  in  air  and  solution  were  thwarted  by  etching  of  these  surfaces.  However,  the 
(001)  surfaces  were  mechanically  robust  toward  tapping  mode  imaging  performed  in 
solution  (Hansma  et  al.,  1994).  This  imaging  mode  measures  changes  in  the  amplitude  of 
an  oscUladng  cantilever  tip,  and  minimizes  the  contact  and  lateral  forces  between  the  tip  and 
the  sample.  Furthermore,  solution  imaging  eliminates  capillary  forces  and,  in  appropriate 
solvents,  reduces  attractive  weaker  van  der  Waals  forces  between  the  tip  and  sample 
(Hutter  and  Bechhoefer,  1993).  The  effective  loading  forces  exerted  during  tapping  mode 
imaging  cannot  be  measured  directly,  but  have  been  estimated  to  be  below  200  pN 
(Radmacher  et  al.,  1994).  These  effects  conspire  to  make  tapping  mode  in  solution  an  ideal 
approach  for  investigating  very  soft  materials,  as  recently  demonstrated  by  observations  of 
isolated  protein  molecules  on  mica  (Fritz  et  al.,  1995). 

The  stability  of  the  (001)  layers  during  tapping  mode  enabled  the  observation  of 

contrast  attributable  to  ordered  insulin  hexamers  in  the  (001)  planes  of  these  polymorphs 
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(Figure  2).  In  the  case  of  the  R6  rhombohedral  polymorph,  the  molecular-scale  real  and 
Fourier  space  AFM  data  acquired  for  the  (001)  layers  conformed  to  a  hexagonal  cell  with  a 
lattice  dimension  of  42.8  ±  6.4  A,  in  agreement  with  the  45.65  A  spacing  of  hexagonally 
packed  nearest  neighbor  insulin  hexamers  along  the  [iTO]  and  [120]  directions,  as 
calculated  from  the  single  crystal  x-ray  structure.  Similarly,  the  molecular-scale  real  and 
Fourier  space  data  acquired  for  (001)  face  of  the  R6  monoclinic  polymorph  clearly  revealed 
a  nearly  square  lattice,  consistent  with  the  symmetry  of  the  (001)  plane  of  this  polymorph. 
Furthermore,  the  dimension  of  this  lattice  was  45.8  A  ±  1.0  A,  in  good  agreement  with  the 
nearest  neighbor  contacts  along  [110]  and  [iTO]  on  the  (001)  face.  The  AFM  unit  cells, 
deduced  from  the  Fourier  spectra,  are  equivalent  alternative  descriptions  of  the  unit  cells 
chosen  for  single  crystal  X-ray  diffraction  structure  solution.  Comparison  of  the 
molecular-scale  contrast  on  (001)  terraces  separated  by  single  steps  revealed  that  the  protein 
layers  were  rotationally  and  translationally  commensurate,  to  within  4°  and  5  A  for  these 
polymorphs  (these  values  are  comparable  to  the  experimental  error).  This  high  degree  of 
commensurism  during  homoepitaxial  growth  is  illustrative  of  good  crystal  quality  with 
respect  to  ordering  of  the  (001)  layers,  which  is  critical  for  diffraction  analysis. 
Furthermore,  the  layer  stacking  deduced  from  data  acquired  on  consecutive  terraces  is 
consistent  with  the  space  groups  of  the  respective  forms,  which  do  not  have  close-packed 
layer  stacking  (Wilson,  1992).  These  data  demonstrate  that  AFM  can  be  used  to  confirm 
space  group  symmetries,  determine  lattice  parameters  and  identify  polymorphs,  thereby 
providing  rudimentary  crystallographic  data  for  crystals  too  small  for  single  crystal 
diffraction  measurements. 

It  is  noteworthy  that  molecular-scale  contrast  has  been  observed  in  contact  mode 

AFM  studies  of  canavalin  and  tobacco  mosaic  virus  in  which  the  crystal  lattice  contains 

layers  of  protein  trimers  with  molecular  weights  of  441  kDa  and  4350  kDa,  respectively 

(Land  et  al.,  1995;  Malkin  et  al.,  1995).  However,  relatively  small  proteins  such  as  insulin 

(hexamer  molecular  weight  of  34.8  kDa)  are  more  likely  to  be  damaged  by  the  large  lateral 
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shear  forces  exerted  by  the  AFM  tip  during  contact  mode  imaging,  suggesting  that  tapping 
mode  imaging  is  the  method  of  choice  for  such  systems.  While  molecular-scale  contrast 
generally  is  not  observed  with  tapping  mode  owing  to  its  poorer  inherent  resolution 
compared  to  contact  mode,  the  achievement  of  molecular-scale  contrast  here  can  be 
attributed  partly  to  the  large  lattice  constants  of  the  (001)  layer. 

[Figure  2] 

Tapping  mode  imaging  in  solution  below  supersaturation  revealed  small  aggregates 
on  the  (001)  layers  of  all  three  polymorphs.  These  aggregates,  illustrated  in  Figure  3  A  for 
the  R6  rhombohedral  polymorph,  were  oriented  along  the  [TlO]  direction  and  were  stable 
indefinitely  during  tapping  mode  imaging  at  concentrations  below  supersaturation.  Their 
morphology  and  shape  orientation  were  independent  of  scan  orientation.  The  aggregates 
appear  to  be  composed  of  multiples  of  poorly  defined  ellipsoids  having  average  lateral 
dimensions  of  110  A  x  52  A  and  a  nominal  height  of  28  A.  This  size,  although  possibly 
inaccurate  due  to  tip  effects,  is  roughly  equivalent  to  2  insulin  hexamers.  Based  on  the 
solution  aggregation  behavior  of  insulin,  it  seems  likely  that  these  features  reflect 
attachment  of  precrystallization  aggregates  formed  in  solution  (Milthorpe  et  al.,  1977; 
Bohinder  and  Geissler,  1984;  Pederson  et  al.,  1994;  Jeffrey,  1974;  Hvidt,  1991).  Similar 
behavior  has  been  reported  for  other  proteins  (Wilson  and  Pusey,  1992;  Pusey,  1991; 
Ataka  and  Asai,  1990;  Durbin  and  Carlson,  1992;  Sazaki  et  al.,  1993;  Kadima  et  al.,  1990) 
and  theoretical  models  of  protein  attachment  kinetics  have  suggested  that  aggregate  addition 
has  at  least  the  same  probability  as  that  for  individual  monomer  units.  Furthermore,  the 
probability  of  multiple  binding  sites  compensates  for  the  decreased  diffusivity  of  the  larger 
aggregates  (Noever,  1995a,  1995b). 

[  Figure  3] 
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If  the  solution  was  allowed  to  reach  supersaturation  by  simply  allowing  a  small 
amount  of  water  to  evaporate,  layer-by-layer  growth  commenced  in  which  (001)  terraces 
grew  by  step  flow  across  the  image  frame  (Figure  4).  Under  these  conditions,  the  terraces 
on  the  insulin  (001)  face  generally  had  areas  exceeding  10  |im2,  with  relatively  few  steps. 
This  suggests  that  the  (001)  face  has  a  relatively  low  surface  energy  and  that  attachment  to 
the  steps  is  thermodynamically  preferred  over  two-dimensional  nucleation.  The  growth 
rate  at  supersaturation  and  room  temperature  was  approximately  10^  nm^  sec^  in  terms  of 
terrace  growth  and  =  5  nm  sec-1  in  terms  of  advancement  of  the  growth  front.  The  growth 
rates  were  calculated  by  comparison  to  the  location  of  stationary  defects  which  served  as 
landmarks.  Careful  examination  of  time  lapse  images  indicates  that  growth  is  fastest  along 
the  [110],  [iTO]  and  [120]  directions  for  the  rhombohedral  forms  and  along  the  [110]  and 
[iTO]  directions  for  the  R6  monoclinic  form.  These  directions  represent  the  closest  contact 
between  insulin  hexamers,  and  accordingly,  would  be  expected  to  be  the  fast  growth 
directions.  From  careful  measurements  of  local  radius  of  curvature  at  multiple  points  along 
the  larger  growing  step  edges  («  131  ±  59  nm),  we  estimate  the  free  energy  of  a  unit  step 
edge,  a  ,  to  be  approximately  3  x  10-'^  J  cm2  with  a  corresponding  kinetic  growth 
coefficient,  p,  to  be  10-3  .  io-4  cm  sec-1  (Burton  et  al.,1951).  These  values  are  comparable 
to  those  determined  from  recent  AFM  studies  of  canavalin  and  tobacco  mosaic  virus  crystal 
growth  (Land  et  al.,  1995;  Malkin  et  al.,  1995),  but  are  approximately  2  orders  of 
magnitude  smaller  than  those  determined  for  inorganic  systems  (Sbhnel,  1982).  The  small 
a  value  is  most  likely  due  to  the  large  solvent  content  of  the  proteins,  which  minimizes  the 
difference  between  the  step  surface  energy  and  the  protein  in  solution.  The  small  p  value 
may  reflect  the  difficulty  in  achieving  the  correct  orientation  of  either  single  proteins  or 
protein  aggregates  during  attachment  at  the  step  sites,  or  the  low  diffusion  coefficient  for 
insulin  (»  7.9  x  lO'^  cm^  sec"!  (Hvidt,  1991))  or  its  precrystallization  aggregates. 
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Under  conditions  of  step  growth,  the  aforementioned  aggregates  typically  are 
consumed  by  the  steps  as  they  flow  across  the  crystal  surface  and  form  new  (001)  layers, 
although  we  cannot  determine  if  the  aggregates  directly  attach  to  the  steps  or  dissociate  to 
hexamers,  dimers  or  monomers  which  subsequently  diffuse  along  the  surface  and  attach  to 
the  step.  Occasionally,  these  aggregates  formed  a  surface  diffusion  barrier  to  growth  at  the 
step  edge,  as  evidenced  by  steps  that  grew  around  the  aggregate,  leaving  a  small  hole 
resembling  a  dislocation  core  which  persists  upon  growth  of  additional  (001)  layers  (Figure 
3B-E).  These  defects  are  similar  to  those  reported  recently  for  lysozyme,  canavalin  and 
tobacco  mosaic  virus  crystals  (Land  et  al.,  1995;  Malkin  et  al.,  1995;  Konnert  et  al.,  1994; 
Durbin  et  al.,  1993;  Durbin  and  Carlson,  1992).  While  these  features  do  not  appear  to 
otherwise  affect  the  growth  of  the  (001)  terraces,  it  is  likely  that  they  occlude  mother 
liquor,  which  can  be  deleterious  for  single  crystal  x-ray  diffraction  and  the  physical 
properties  of  commercial  formulations.  These  defects  differ  from  pits  introduced  by 
intentional  mechanical  etching  with  the  AFM  tip,  which  are  transient  and  rapidly  fill  in  by 
the  same  layer  growth  mechanism  described  above.  These  observations  suggest  that  some 
insulin  aggregates  attached  to  the  (001)  face  are  strongly  bound  and  cannot  achieve  the 
correct  orientation,  migrate  to  the  step,  or  dissociate  into  more  mobile  species  on  the  time 
scale  of  the  terrace  growth,  thereby  instigating  the  creation  of  dislocation  cores. 

Previous  contact  mode  AFM  studies  of  crystal  growth  have  indicated  that  apparent 

growth  rates  are  faster  in  the  region  being  scanned  compared  to  surrounding  regions.  This 

may  be  due  to  defect  formation  caused  by  the  tip  (LaGraff  and  Gewirth,  1995,  1994)  or 

enhanced  convection  caused  by  the  rastering  motion  of  the  tip.  After  extensive  searching, 

we  have  not  found  any  evidence  for  tip-induced  defect  formation  under  typical  scanning 

conditions  nor  have  we  detected  any  screw  dislocations  that  may  serve  as  the  source  of 

steps,  suggesting  very  low  densities  of  these  growth  centers.  Rather,  growing  terraces 

entered  the  scanning  frame  from  outer  regions,  indicating  that  the  tip  was  not  initiating 

growth.  After  the  growing  terrace  filled  the  scanning  frame,  expanding  the  field  of  view 
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about  the  original  scanned  region  revealed  that  growth  was  somewhat  faster  in  the  scanned 
regions  than  in  the  oudying  regions,  similar  to  behavior  reported  for  canavalin  (Land  et  al., 
1995).  These  observations  suggest  that  a  very  small  number  of  dislocations  initiate  the 
formation  of  steps  remote  from  the  scanning  region,  but  that  the  tip  enhances  mixing  of  the 
boundary  layer  at  the  growing  crystal  surface  as  a  consequence  of  its  vertical  and  horizontal 
motion,  thereby  enhancing  otherwise  diffusion  limited  growth. 

[Figure  4] 

These  studies  illustrate  that  tapping  mode  AFM  performed  directly  in  solution 

provides  resolution  of  molecular-scale  contrast  which  enables  identification  of  different 

polymorphs  of  bovine  insulin,  confirmation  of  the  crystal  packing  and  symmetry,  and 

insight  into  the  growth  mechanisms.  The  apparent  absence  of  screw  dislocations,  although 

the  actively  growing  steps  most  likely  emanate  from  them,  indicate  that  insulin  crystals  can 

be  grown  with  a  minimal  number  of  these  defects.  However,  impurities  or  misaligned 

aggregates  bound  to  actively  growing  crystal  faces  can  instigate  the  formation  of  core 

dislocations.  We  anticipate  that  the  ability  to  resolve  molecular  scale  contrast  while 

simultaneously  visualizing  the  nanoscale  topography  and  growth  modes  of  the 

morphologically  important  faces  of  insulin  single  crystals  will  become  a  useful  approach 

for  optimizing  crystallization  conditions  of  this  therapeutically  important  protein  and  its 

modified  forms.  The  capabilities  of  in  situ  AFM  provide  an  opportunity  for  direct 

examination  of  the  influence  of  various  experimental  parameters  during  growth  at  the  near- 

molecular  level,  instead  of  relying  on  inference  from  inspection  of  mature  crystals.  This 

includes  studies  aimed  toward  elucidation  of  the  role  of  additives  on  polymorph  selectivity 

(Kaarsholm  et  al.,  1989;  Derewenda  et  al.,  1989,  Smith  and  Dodson,  1992,  Ciszak  and 

Smith,  1994),  solid  state  structure  (Smith  and  Ciszak,  1994)  and  crystal  morphology 

(Markman  et  al.,  1992)  protein-receptor  interactions  (Lee  et  al.,  1994a,  1994b  ;  Boland 
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and  Ratner,  1995;  Florin  et  al.,  1994;  Lee  et  aL,  1994;  Moy  et  al.,  1994;  Chilkoti  et  al., 
1995;  Ohnseorge  et  al.,  1992;  Drake  et  al.,  1989;  Weisenhom  et  al.,  1992),  and  the 
formation  of  cocrystals  with  other  proteins  (Low  et  al.,  1974;  Krayenbiihl  and  Rosenberg, 
1946;  Hagedom  et  al.,  1936). 
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Figure  Captions 

Figure  1.  (A)  Contact  mode  AFM  image  revealing  large  molecularly  flat  terraces 
observed  on  the  (001)  face  of  R6  rhombohedral  bovine  insulin  in  mother  liquor.  Terrace 
heights  are  consistent  with  the  c-axis  lattice  constant  for  this  polymorph  (c  =  40.39  A). 
Scan  size  :  2.5  itm  x  2.5  itm.  Scan  rate  :  5.55  Hz.  (B)  Tapping  mode  AFM  image  of  the 
(001)  face  of  R6  rhombohedral  bovine  insulin  acquired  in  air,  revealing  large  (==  60  nm) 
rounded  aggregates  on  the  terraces  and  step  edges.  Narrow  canyons  between  adjacent 
terraces  support  a  small  amount  of  lateral  contraction  during  fracturing  that  accompanies 
crystal  dehydration.  Terrace  heights  obtained  in  air  are  consistent  with  the  thickness  of  an 
isolated  insulin  hexamer  unit.  Scan  size  :  2.5  |im  x  2.5  ^im.  Scan  rate  :  4.60  Hz.  Inset: 
Terrace  and  ledges  of  T6  rhombohedral  polymorph  imaged  in  air  under  tapping  mode 
conditions.  The  steps  are  scalloped  and  have  well  defined  120°  angles,  consistent  with 
(hkO)  step  planes,  most  likely  (110).  Inset  scan  size  :  1.25  |J.m  x  1.25  |im.  Scan  rate  :  12 

Hz.  Height  profiles  are  taken  along  lines  marked  A  -  A' 

Figure  2.  (A)  Raw  and  Fourier  filtered  real  space  data  for  the  (001)  face  of  the  R6 
rhombohedral  polymorph  of  bovine  insulin.  Scan  size  :  50  nm  x  50  nm.  Scan  rate  :  5.55 
Hz.  Inset:  Two-dimensional  Fourier  power  spectrum  revealing  reciprocal  lattice  vectors 
corresponding  to  real  space  periodicities  of  ai  =  41.5  A;  32  =  37.2  A;  33  =  49.8  A;  712 
=  47°;  723=  54°-  The  slight  distortion  from  hexagonal  symmetry  is  attributed  to  scanning 
artifacts.  (B)  Raw  and  Fourier  filtered  real  space  data  for  the  (001)  face  of  the  R6 
monoclinic  polymorph  of  bovine  insulin.  Inset:  Two-dimensional  Fourier  power  spectrum 
revealing  reciprocal  lattice  vectors  corresponding  to  real  space  periodicities  of  3i  =  46.4 
A;  32  =  45.1  A;  Yi2  =  84°.  Scan  size  :  50  nm  x  50  nm.  Scan  rate  :  8.14  Hz.  (C)  Insulin 
hexamer  packing  motif  of  the  (001)  face  of  the  R6  rhombohedral  polymorph  of  bovine 
insulin  based  on  single  crystal  X-ray  diffraction  data.  The  real  space  vectors  3i ,  32  and 
33  correspond  to  the  nearest  neighbor  hexamer  vectors  [120],  [210],  [iTO],  respectively. 
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Figure  Captions  (cont'd) 

Figure  2.  (cont'd) 

(D)  Insulin  hexamer  packing  motif  of  the  (001)  face  of  the  R6  monoclinic  polymorph  of 
bovine  insulin  based  on  single  crystal  X-ray  diffraction  data.  The  real  space  vectors  ai 
and  32  correspond  to  the  nearest  neighbor  hexamer  vectors  [110]  and  [iTO],  respectively. 
Unit  cells  based  on  single  crystal  X-ray  diffraction  data  are  shown  as  solid  black  lines. 
The  R6  rhombohedral  polymorph  crystallizes  in  space  group  R3  with  lattice  parameters  a  = 
b  =  79.92  A,  c  =  40.39  A  indexed  on  hexagonal  axes.  The  R6  monoclinic  polymorph 
crystallizes  in  space  group  P2i  with  lattice  parameters  o  =  61.23  A,  Z?  =  61.65  A,  c  =  48.05 
A,  P  =  110.50°.  AFM  unit  cells  corresponding  to  the  two-dimensional  Fourier  power 
spectra  are  depicted  as  solid  white  lines.  Note  that  the  position  of  the  AFM  unit  cells  in  the 
x-y  plane  is  arbitrary.  The  AFM  unit  cells,  deduced  from  the  Fourier  spectra,  are 
equivalent  alternative  descriptions  of  the  unit  cells  chosen  for  single  crystal  X-ray 
diffraction  structure  solution.  Structural  models  prepared  using  Cerius^  version  1.6 
(BIOSYM  /  Molecular  Simulations  Inc.,  San  Diego,  CA). 

Figure  3.  (A)  In  situ  AFM  tapping  mode  image  of  protein  aggregation  on  the  (001)  face 
of  the  R6  rhombohedral  polymorph  of  bovine  insulin.  The  aggregates  appear  to  be 
oriented  along  the  [iTO]  direction,  as  deduced  from  comparison  to  the  molecular-scale 
contrast  of  the  underlying  (001)  crystal  face.  Scan  size  :  300  nm  x  3(X)  nm.  Scan  rate  :  2.97 
Hz.  (B)-(E)  Real-time  in  situ  tapping  mode  images  collected  at  0,  643,  772  and  801 
seconds,  respectively,  revealing  formation  of  persistent  core  dislocations  instigated  by 
misaligned  protein  aggregrates.  Note  that  one  complete  terrace  has  filled  the  scanning  frame 
between  images  B  and  C.  A  new  terrace  has  entered  the  frame  in  image  C  and  has  filled 
the  scanning  frame  in  image  E.  Outlines  indicate  location  of  aggregates  and  core 
dislocations.  Image  size:  1.5  |J,m  x  3.5  pm.  Scan  rate  :  2.0  Hz. 
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Figure  Captions  (cont'd) 

Figure  4.  (A)-(D)  Real-time  in  situ  tapping  mode  images  of  R6  rhombohedral 
polymorph  of  bovine  insulin  collected  at  0,  257,  644  and  1287  seconds,  respectively, 
revealing  layer-by-layer  growth  of  individual  insulin  terraces.  Scan  size:  5  |im  x  5  |im. 
Scan  rate  :  2.0  Hz.  The  arrow  indicates  the  position  of  a  stationary  defect  site  used  as  a 
landmark  for  calculating  terrace  growth  rates. 
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